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The Light-front quark model (LFQM) has been applied to calculate the transition matrix elements 
of heavy hadron decays. However, it is noted that using the traditional wave functions of the LFQM 
{3J[), given in literature, the theoretically determined decay constants of the T{nS) obviously contradict 

^ . to the data. It implies that the wave functions must be modified. Keeping the orthogonahty 

' among the nS states and fitting their decay constants we obtain a series of the wave functions 

I for T(nS'). Based on these wave functions and by analogy to the hydrogen atom, we suggest a 

■ modified analytical form for the T{nS) wave functions. By use of the modified wave functions, 

the obtained decay constants are close to the experimental data. Then we calculate the rates of 
radiative decays of T{nS) — >■ r;;, -I- 7. Our predictions are consistent with the experimental data on 
decays T{3S) ^ 771, + 7 within the theoretical and experimental errors. 



Q^' PACS numbers; 13.25.Gv, 13.30.Ce, 12.39.Ki 



I. INTRODUCTION 



0^ ■ Since the relativistic and higher-order as corrections are less important for bottomonia than for any other qq 
■ systems, study on bottomonia may offer more direct information about the hadron configuration and apphcation of 
the perturbative QCD. The key problem is how to deal with the hadronic transition matrix elements which are fully 
\G • governed by the non-perturbative QCD effects. Many phenomenological models have been constructed and applied. 
[ Each of them has achieved relative successes, but since none of them are based on any well established underlying 
theories, their model parameters must be obtained by fitting data. By doing so, some drawbacks of the model are 
exposed when applying to deal with different phenomenological processes. Thus one needs to continuously modify the 
model or re-fit its parameter, if not completely negate it. The light front quark model is one of such models. It has 
been applied to calculate the hadronic transitions and generally considered as a successful one. The model contains 
^ I . a Gaussian-type wavefunction whose parameters should be determined in a certain way. 
d ' The Gaussian- type wavefunction was recommended by the authors of Refs. [l|, 12| and most frequently the wave- 
function for harmonic oscillator is adopted which we refer as the traditional LFQM wavefunction. As we employed the 
traditional LFQM wave functions to calculate the branching ratios of T(nS) — > 77b -I- 7, some obvious contradictions 
between the theoretical predictions and experimental data emerged. Namely, the predicted B{T{2S) ?/h + 7) was 
one order larger than the experimental upper bound 0]. Moreover, as one carefully investigates the wave functions, 
he would face a serious problem. If the traditional wave functions were employed, the decay constants of T(nS) (/y) 
would increase for higher n. It obviously contradicts to the experimental data and the physics picture which tells us 
that the decay constant of a nS state is proportional to its wavefunction at origin which manifests the probability that 
the two constituents spatially merge, so for excited states the probability should decrease. Thus the decay constants 
should be smaller as n is larger. The experimental data confirm this trend. But the theoretical calculations with the 
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traditional wave functions result in an inverse order. To overcome these problems, one may adopt different model 
parameters (refers to (3) by fitting individual n's decay constants as done in 0, 1^, but the orthogonality among 
the nS states is broken. In this work, we try to modify the harmonic oscillator functions and introduce an explicit 
n-dependent form for the wave functions. Keeping the orthogonality among the nS states (n = 1, ...5), we modify the 
LFQM wave functions. By fitting the decay constants of T{nS), the concerned model parameters are fixed. Besides 
fitting the decay constants of the T(nS) family, one should test the applicability of the model in other processes. We 
choose the radiative decays of T{nS) — >■ Tyf, + 7 as the probe. As a matter of fact, those radiative decays are of great 
significance for understanding the hadronic structure of bottomonia family. 

Indeed, the spin-triplet state of bottomonia T{nS) and the P-states Xb{nP) were discovered decades ago, however 
the singlet state rj^ evaded detection for a long time, even though much efforts were made. Many phenomenological 
researches on rji, have been done by some groups [5l-[l2|. Different theoretical approaches result in different level 
splitting AM = T (15) - 77b (IS*). In Q the authors used an improved perturbative QCD approach to get AM = 44 
MeV; using the potential model suggested in Eichten and Quigg estimated AM = 87 MeV @; in Ref. 0] 
the authors selected a non-relativistic Hamiltonian with spin dependent corrections to study the spectra of heavy 
quarkonia and got AM =57 MeV; the lattice prediction is AAf=51 MeV Q, whereas the lattice result calculated in 
Ref. Q was AM = 64± 14McV. Ebcrt etal. directly studied spectra of heavy quarkonia in the rclativistic quark 
model and gave m^j_ = 9.400 GeV. The dispersion of the values may imply that there exist some ambiguities in our 
understanding about the structures of the bb family. 
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FIG. 1: AM coming from different experimental measurement and theoretical work. 



The Babar Collaboration Q first measured B{T{3S) -frit) = (4.8 ± 0.5 ± 0.6) x 10"'', and determined m,,^ = 
9388.9t^;3 ± 2.7 MeV, AM = 71.41^ ;^ ± 2.7 MeV in 2008. New data m^, = 9394.2t^J ± 2.0 MeV and B{T{2S) 
rjrjt) = (3.9 ± l.ltll) X 10"'^ were released in 2009 [H. More recently the CLEO Collaboration ^ confirmed the 
observation of rji, using the database of 6 million T(3S') decays and assuming r(77f,) «10 MeV, they obtained B{T{iS) — > 
^nb) = (7.1 ± 1.8 ± 1.1) X 10"'', m^, = 9391.8 ± 6.6 ± 2.0 MeV and the hyperfine splitting AM = 68.5 ± 6.6 ± 2.0 
MeV, whereas using the database with 9 million T(25) decays they obtained ;B(T(2S') — > "fq^) < 8.4 x 10^^ at 90% 
confidential level. It is noted that the data of the two collaborations are in accordance on m^j, , but the central values 
of B{T{3S) — ^rjb) arc different. However, if the experimental errors are taken into account, the difference is within 
one standard deviation. 

Some theoretical works [l7l - [T9l | are devoted to account the experimental results. In Ref. [13] the authors studied 
these radiative decays and estimated B{T{3S) r]t + -f) = 4 x 10"", 6(T(25) ^ rjt + -f) = 1.5 x lO"'^ and 
S(T(1S') -> 77b + 7) = 1.1 X 10"'' with the mass m^, = 9.400 GeV. Their results about to^, and S(T(3S') Vb+l) are 
close to the data. The authors of Ref. ^20^ systematically investigated the magnetic dipole transition V — >■ P7 in the 
light-front quark model (LFQM) jX, 2, 21, 22]. In the QCD-motivated approach there are several free parameters, i.e., 
the quark mass and /3 in the wave function (the notation of /3 was given in the aforementioned literatures) which are 
fixed by the variational principle, then B{T{1S) rji, + j) was calculated and the central value is 8.4 (or 7.7) x 10~^ 
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^. It is also noted that the mass of m^^ = 9.657 (or 9.295) GeV presented in Ref. [1^ deviates from the data hsted 
before, so we are going to re-fix the parameter /3 in other ways namely we fix the parameter (3 by fitting data. 

Since experimentally, m^j, is determined by B{T{nS) — > + 7) and a study on the radiative decays can offer us 
much information about the characteristics of jyt,, one should carefully investigate the transition within a relatively 
reliable theoretical framework. That is the aim of the present work, namely we will evaluate the hadronic matrix 
element in terms of our modified LFQM. 

This paper is organized as follows. After this introduction, in section II we discuss how to modify the traditional 
wave functions in LFQM. We present the formula to calculate the form factors for V — !• Pj in the LFQM and 
numerical results in section III. The section IV is devoted to our conclusion and discussion. 



II. THE MODIFIED WAVE FUNCTIONS FOR THE RADIALLY EXCITED STATES 



When the LFQM is employed to calculate the decay constants and form factors, one needs the wave functions of the 
concerned hadrons. In most cases, the wave functions of harmonic oscillator are adopted. In the works [H. [2|. [20l - [23j . 
only the wave function of the radially ground state is needed, but when in the processes under consideration radially 
excited states of are involved, their wave functions should also be available. In [23, HI], the traditional wave functions 
if for IS* and 25* states in configuration space from harmonic oscillator are given as 



./32 3/4 



.-(0 ^ (^)-exp(4,V)-L(3-2,V) 



(1) 



In order to maintain the orthogonality among nS states, the parameter /? in the above two functions are the same. 
The wave functions for other nS state can be found in Appendix A. 

The decay constants of the nS states are directly proportional to the wave function at the origin 



fv oc ip{r = 0). 



(2) 



If we simply adopt the wave functions of harmonic oscillator for all of them as we do for the 15 state, then we find the 
wave functions at the origin, i.e. (p{r = 0) (see Appendix for details) rises with increase of n (the principle quantum 
number) which means the decay constants would increase for larger n. For example, by Eq. ([!]) the ratio of wave 
functions of 25" and IS states at the origin is 3/-\/6 > 1. 

The decay constants fv of T{nS) are extracted from the processes r(T(nS') — > e+e^) with 
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(3) 



where V represents T{nS) and my its mass. By use of the experimental data from PDG [26| . we obtain the 
experimental values for /y which are listed in Table ID Obviously, the decay constant becomes smaller as n is larger. 
In LFQM, the formula for calculating the vector meson decay constant is given by [l], Q 



fv 
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dx I d k { 



(t>{nS) 



^2x{l- x)Mo 



xAU 



mi(mi — m2) — k\ + 



mi + m2 
/ 

Mq + mi + 1712 



(4) 



where mi = m2 = ?7it and other notations are collected in the Appendix. In the calculation we set m^ = 5.2 GeV 
following [2^ and the decay constant of T(IS') is used to determine the parameter /3t as the input. Wc obtain 
/3t = 1.257±0.006 GeV corresponding to f^^g^ = 715 ± 5 McV. In order to illustrate the dependence of our results 
on mf), we re-set m^ = 4.8 GeV to repeat our calculation, then by fitting the same data, we fix /3x = 1.288±0.006 GeV 
and all the results are clearly shown in the following tables. The in Table U are the decay constant calculated in 
the traditional wave functions. These results expose an explicit contradictory trend. Thus, our calculation indicates 
that if the traditional wave functions arc used, the obtained decay constants of T{nS) would sharply contradict to 
the experimental data. 



^ The difTorcnt values correspond to the different potentials adopted in the calculations. 
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TABLE I: The decay constants of T(7iS) (in the unit of MeV). The column "/x" represents the theoretcal predictions with the 
traditional wave function in LFQM. The column "/^" represents the prediction with our modified wave function and the values 
in the brackets are the corresponding values with nib = 4.8 GeV as input. (The other values are corresponding to mt — 5.2 
GeV.) 



nS 


rcxp 
Jt 


/I 


Jr 


IS 


715±5 


715±5 


715±5 (715±5) 


2S 


497±5 


841±7 


497±5 (498±5) 


3S 


430±4 


925 ±8 


418±5 (419±4) 


4S 


340±19 


993 ±8 


378±4 (397±4) 


5S 


369±42 


1040 ±9 


349±4 (351±4) 



As aforementioned, the wave functions must be modified. Our strategy is to establish a new Gaussian-type wave 
function which is different from that of harmonic oscillator. As modifying the wave functions, several principles must 
be respected: 

(1) The wave function of IS" should not change because its application for dealing with various processes has been 
tested and the results indicate that it works well; 

(2) The number of nodes of nS should not be changed; 

(3) A factor may be added into the wave functions which should uniquely depend on n in analog to the wave 

function of the hydrogen-like atoms which is written as Rn{r) = P^'^'^'^{r)e~~ , where Pnir) is a polynomial and Z 
is the atomic number, ao is the Bhor radius; 

(4) Using the new Gaussian-type wave function, the contradiction for the decay constants can be solved. 

In the LFQM, we only need the wave functions in the momentum space. Fourier transformation gives us the 
corresponding forms in the momentum space, see the Appendix for details. The 15 wave function is remained and 
used to fix the model parameter. Now let us investigate the wave function of 25*. According to the analog to the 
hydrogen-like atom, we introduce a factor g2 represents n-dependence to the exponential in the wave function of 25, 
thus the wave function of 2S is changed to 




where the subscript M denotes the modified function. Then by requiring it to be orthogonal to that of 15 and 
normalizing the wave function, we determine the parameters a and b in the modified wave function of 25. With this 
new wave function of 25, we demand the theoretical decay constant be consistent with data so 52 should fall into a 
range determined by the experimental errors. Going on, we obtain the modified wave function of 35 and that for 4S 
and 5S as well. In this case the modified wave functions of nS states are more complicated than the traditional ones. 

We have gained a series of numerical gnS by the principles we discussed above, then we wish to guess an analytical 
factor gn which is close to the numerical values of the series. We find that if g„ = S = 1/1.82) is set, we almost 
recover the numerical series. Thus the wave function of the nS state in the momentum space can be written as 

Cf(p') = Pn(p')cxp(-n*^), (6) 

where P„(p^) is a polynomial in p^. The corresponding wave function of the nS state in the configuration space can 
be written as 

rJ{r)=P:Ar')e^p{-^). (7) 
Comparing with the case of the hydrogen-like atoms which the n5-wavc functions arc written as 

i?„oM = P„^-^'^'^(r)exp(^). (8) 

in the configuration space, where P^^'^''(r) is a polynomial in r. The factor \ jn in the exponential power is obtained 
by solving the Schrodinger equation where only the Coulomb potential exists. To modify the wave functions we 
get the factors numerically for all the nS states, then "guess" its analytical form. In the LFQM, the factor Xjn^ is 
introduced to fit the experimental data for nS decay constants. Definitely this analytical form is not derived from 
an underlying theory, such as that for the hydrogen atom, thus the dependence on n is only an empirical expression. 
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But we are sure that if the model is correct and our guess is reasonable, it should be obtained from QCD (maybe 
non-perturbative QCD). It is noted that the experimental errors are large, so that other forms for g„ might also be 
possible. The theoretical estimation of the decay constants of T{nS) (fj^) are also presented in Table HI The modified 
wave functions seem to work well and they could be used for evaluating B{T(nS) ^ 776 + 7). 



III. THE TRANSITION OF T{nS) Vb + 1 



In this section, we calculate the branching ratios of T{nS) rib + ^ in terms of the modified wave functions derived 
in the above section. 



A. Formulation of TinS) ^ + 7 in the LFQM 

The Feynman diagrams describing T{nS) — > 776+7 are plotted in Fig. [51 The transition amplitude of T(n5') — > 77b +7 
can be expressed in terms of the form factor J^T{nS)-¥rn,{l^) which is defined as (20I [2H 

{rib{V')\JliJT{V.h)) ^ 7ee^'''"^e.(P,/7)gp7'.J-T(„5)-.rJg'), (9) 

where V and V' are the four-momenta of T{nS) and 77;,. q ^ V is the four-momentum of the emitted photon and 
e^{V, h) denotes the polarization vector of T{nS) with helicity h. For applying the LFQM, we first let the photon be 
virtual, i.e. leave its mass-shell ~ into the un-physical region of q^ < 0. Then ^T(riS)-!-?7b('?^) can be obtained in 
the (7+ = frame with q^ = q^q^ ^ = < 0- Then we just analytically extrapolate -7^T(riS)-^rft (qI) from the 
space-like region to the time-like region (q^ > 0). By taking the limit q"^ — >■ 0, one obtains J-T{nS)~>nbi^'^ — 




w U w , -^''^ W 

FIG. 2: Feynman diagrams depicting the radiative decay T(nS') — >■ ryt + 7. 



By means of the light front quark model, one can obtain the expression of the form factor J^T{nS)^7^bil^) po} - 

^T(nS)^nM^) = ebl{mi,m2,q'^) + ebl{m2,mi,q'^), (10) 
where Cb is the electrical charge for the bottom quark, mi = 777,2 = 771;, and 



^ dx 



{x,'k' ±)(j){x,'k±) 
xiMoM' 



' Mo qi 



(11) 



where A = X2mi + Xim2, x = xi and the other variables in Eq. pip are defined in Appendix. In the covariant 
light-front quark model, the authors of [2l| obtained the same form factor J^x(nS)^?)i, (q^)- The decay width for 
T(7iS') —> 776 -I- 7 is easily achieved 



r(T(7i5)^7,b+7)-- 



2771 



T{nS) 



J"- 



T(nS)->r,i, 



(0). 



(12) 



where a is the fine-structure constant and 771x(tiS): "^jyb ''^^^ the masses of T{nS) and rjb respectively. 
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TABLE II: The branching ratios of T{nS) 77)6. In the column "Sf ", mt = 5.2GeV, /3t = 1.257 ± 0.006 GeV and 
I3nt = 1.246 ± 0.005 GeV. In the column "Bff", mt = 4.8GeV, /3t = 1-288 ± 0.006 GeV and Z?^, = 1.287 ± 0.005 GeV. In the 
column "B'^", nib = 5.2GeV, /3t = 1.257 ± 0.006 GeV and /3^, = 1.249 ± 0.005 GeV. 



Decay mode 














Experiment 




T(IS') - 




(1.94 ± 0.41) X 10" 


-4 


(2.24 ±0.47) X 10" 


-4 


(1.94 ±0.42) X 10" 


-4 






T(2S') - 


-)• 77b + 7 


(3.90 ± 1.49) X 10" 


-4 


(3.90 ± 1.49) X 10" 


-4 


(3.90 ± 1.49) X 10" 


-4 


(3.9± l.l^o g) X 10- 


' [15J 


T(3S) - 


776 + 7 


(1.87 ±0.71) X 10" 


-4 


(1.68 ±0.72) X 10" 


-4 


(1.05 ±0.40) X 10" 


-b 


(4.8 ±0.5 ±0.6) X 10 
(7.1 ± 1.8 ± 1.1) X 10 


-* [16] 


T(4S) - 


771, + 7 


(8.81 ±3.32) X 10" 


-8 


(7.82 ± 3.35) X 10" 


-8 


(2.25 ± 0.88) X 10" 


-iU 






T(5S') - 


776 +7 


(1.17 ±0.43) X 10" 


-8 


(1.02 ±0.45) X 10" 


-8 


(1.57 ±0.52) X 10" 


-L'2 







B. Numerical results 

Now we begin to evaluate the transition rates of T(2S') — > 77c + 7 with the modified wave functions. We stiU use 
the values of mb = 5.2 GeV and /3x = 1.257 ± 0.006 GeV given in last section. The parameter /3^^ is unknown, we 
determine it from T(2S') — ^ ■yrji, process. Comparing with the data B{T{2S) — > 777;,) = 3.9 x lO""* [lH, we obtain 
/3jj^ = 1.246 ± 0.005 GeV which is consistent with our expectation, namely it is close to the value of (3r = 1.257 GeV. 
Under the heavy quark limit, they should be exactly equal, and the deviation must be of order 0(l/mb) which is 
small jl^. With these parameters, we can calculate the branching ratios B{T{1S) -> i]b + 7), B{T{3S) — > 776 ± 7), 
S(T(4S') — >■ 77(, + 7), and B{T{5S) — >■ 77^ + 7). The numerical results are presented in the column "^Sj^" of Table 
im Indeed, the b-quark mass is an uncertain parameter which cannot be directly measured and in some literatures, 
different values for b-quark mass have been adopted. To see how sensitive to the b-quark mass the result would be, 
we also present the numerical resuhs with mt = 4.8 GeV, /3t = 1.288 ± 0.006 GeV and P^, = 1.287 ± 0.005 GeV 
in the column "B^" of Table HIl The results in the column "6"^" of Table HIl are obtained with the traditional wave 
functions. Apparently, as the modified wave functions are employed, the theoretical predictions on branching ratios of 
the radiative decays are much improved, namely deviations from the data is diminished. About the numerical results, 
some comments are given as following: 

(1) Comparing the results shown in column Bf^ with those in column B^, we can find that they are not sensitive 
to mb- 

(2) For the decay T(IS') — > 776 + 7, our prediction on the branching ratio is about 2.0 x 10~^. This mode should be 
observed soon in the coming experiment. Our prediction is consistent with the results of Ref. [20j. The branching 
ratio is not sensitive to /3^(,, but sensitive to the mass splitting AM. That is easy to understand. Since the decay 
width is proportional to (AM)^, thus as AM is small, i.e., the masses of initial and daughter mesons are close to 
each other, any small change of ?7i^j_ can lead to a remarkable difference in the theoretical prediction on the branching 
ratio. Thus the accurate measurement on B{T{1S) ^ 776 ± 7) wiU be a great help to determine the mass of m^^. 

(3) The process of T{2S) — >■ ?7f, + 7 is used as an input to determine the parameter of rjb- The prediction of 
T(35') -± r/f, +7 is in accordance with the experimental data by the order of magnitude. After taking into account the 
experimental and theoretical errors, they can be consistent. This result could be of relatively large errors, because 
we only use four parameters {mb, /3t, /^r/j, ct) to determine five decay constants and three branching ratios for 
T{1S, 2S, 3S) r]b + J and ah of them possess certain errors. 

(4) The branching ratios for the processes T(4S') i]b + j and T{5S) rjb+'-f are at the order of 10"®, it is nearly 
impossible to be observed in the near future if there aren't other mechanisms to enhance them. 

(5) As an application, we predict the decay constant of rjb in terms of the model parameters we obtained above. We 
calculate the branching ratio of B{T{2S) — > 77;;,) in the LFQM. By fitting data we fix the concerned model parameters 
for rjb, and then with them we predict the decay constant of rjb in the same framework of the LFQM P, [|. In the 
calculations, b-quark mass mb and are input parameters. 

To show how sensitive the results are to the parameters, we use the two sets of input parameters given above, 
and the corresponding results are as follows: f,j^ = 567 MeV when mb = 5.2GeV and = 1.246GeV ; fjj^ = 604 
MeV when mb = 4.8GeV and = 1.287GeV. For a comparison, we deliberately change only mb while keeping 
unchanged to repeat the calculation and obtain f,f^ — 591 MeV when 7?T,f, = 5.2GeV and = 1.287GeV. It is noted 
that frj^ is more sensitive to /3^^ , rather than mb ■ 
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IV. CONCLUSION 

The LFQM has been successful in phenomenological appUcations. It is beUeved that it could be a reasonable model 
for dealing with the hadronic transitions where the non-perturbative QCD effects dominate. However, it seems that 
the wave function adopted in the previous literature has to be modified. As we study the decay constant of T{nS), 
we find that there exists a sharp contradiction between the theoretical prediction and data as long as the traditional 
harmonic oscillator wave functions were employed. Namely, the larger n is, the larger the predicted decay constant 
would be. It is obviously contradict to the physics picture that for higher radially excited states, the wave function 
at origin should be smaller than the lower ones. But the old wave functions would result in an inverse tendency. If 
enforcing all the decay constants of T{nS) to be fitted to the data in terms of the traditional wave functions, the 
orthogonality among all the nS states must be abandoned, but it is not acceptable according to the basic principle of 
quantum mechanics. 

Thus we modify the wave functions of the radial excited states based on the common principles. Namely, we keep 
the orthogonality among the wave functions and their proper normalization. Moreover, we require the wave functions 
fMij") at origin r = to be consistent with the data, i.e. the decay constants for higher n must be smaller than that 
of the lower states. Concretely, we modify the exponential function in the wave functions by demanding the power 
not to universal for all n's, but be dependent on n. Concretely we add a numerical factor gn into cxp{gn^^) and by 
fitting the data of the decay constants of T{nS) we obtain a series of the numbers of gn- Within a reasonable error 
range, we approximateg„ as gin) = and calculate the value for S. It is an alternative way which is different from 
that adopted in Ref. [20|, to fix the parameter. 

With the modified wave functions of T(nS), we calculate the branching ratios of T{nS) ^ r/b + 7 in the LFQM. 
First by fitting the well-measured central value of B{T{2S) — >■ ?7h + 7) [lal; we obtain the parameter By the 
effective heavy quark theory, in heavy quark limit the spin singlet and triplet of bb system should degenerate, namely 
the parameters of /3t(is) and /^r/b should be very close. Our numerical result confirms this requirement. 

Then we estimate the other T(nS) — > jyti + 7. The order of magnitudes of our numerical results is consistent with 
data. Even though the predicted branching ratios still do not precisely coincide with the data, the result is much 
improved. The branching ratios of processes T{AS) — >■ f?6 + 7 and T{5S) — >■ ??& + 7 are predicted to be at the order 
of 10~^. They are difficult to be measured in the future as long as there is no new physical mechanism to greatly 
enhance them. 

By studying the radiative decay of T(nS) — > + 7, we can learn much about the hadronic structure of r]b. Even 
though much effort has been made to explore the spin singlet rjt, in particle data group (PDG) of 2008, rji, was still 
omitted from the summary table j26| . In fact, determination of the mass of rjb is made via the radiative decays of 
T{nS) r/b + 7 [13, and the recent data show = 9388. 9t 2 l{stat) ± 2.7{syst) MeV by the T{3S) data and 
rrin^ = 9394.2;4 g(stai) ± 2.0{syst) MeV by the T(2S') data [H. Penin [H reviewed the progress for determining 
the mass of r]b and indicated that the accurate theoretical prediction of m,,j, would be a great challenge. Indeed, 
determining the wave function of rjb would be even more challenging. We carefully study the transition rates of the 
radiative decays which would help to extract information about Wr/b. The transition rate of T(IS') — >■ rjb + 1 is very 
sensitive to the mass splitting AM = ?nT(iS) ~ ™r;!, due to the phase space constraint, thus an accurate measurement 
of the radiative decay may be more useful to learn the spin dependence of the bottominia. 
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Appendix 



A. The radial wave functions 



The traditional wave functions (j) in configuration space from harmonic oscillator [2J| are 
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and their Fourier transformation are 
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The modified wave functions tpM in configuration space are defined 
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with coefficients, which are irrational numbers and are kept five digits after the decimal point 
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The corresponding modified wave functions in momentum space are 



with, coefficients 
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B. Some notations in LFQM 

The incoming (outgoing) meson in Fig. [2]has the momentum P^'^ = pi^'' +P2 where pi^''' and p2 are the momenta 
of the off-shell quark and antiquark and 
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with xi + X2 — 1, where Xi and kj_(k'^) are internal variables. Mq and Mq are defined 
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More information can be found in Rcf. 
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